Disorder by order in graphene 
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We predict the existence of an intriguing "disorder by order" phenomenon in graphene transport 
where higher-quahty (and thus more ordered) samples, while having higher mobility at high carrier 
density, wiU manifest more strongly insulating (and thus effectively more disordered) behavior as 
the carrier density is lowered compared with lower quality samples (with higher disorder) which 
exhibit an approximate resistivity saturation phenomenon at low carrier density near the Dirac 
point. This predicted behavior simulating a metal-insulator transition, which we believe to have 
recently been observed in an experiment at Manchester University^, arises from the suppression of 
Coulomb disorder induced inhomogeneous puddles near the charge neutrality point in high quality 
graphene samples. 



I. INTRODUCTION 

An electronic material, metal or doped semiconductor, 
typically exhibits higher low-temperature conductivity as 
the amount of quenched random disorder is decreased in 
the system, i.e., as the system becomes more ordered. 
It is therefore a universal expectation that a purer metal 
with lower impurity disorder would always exhibit higher 
low-temperature conductivity than a dirtier metal with 
higher disorder. 

In the current work, we theoretically establish a 
counter-intuitive possibility in graphene which is in sharp 
contrast to the universal scenario of increasing conduc- 
tivity with increasing order. We show that in monolayer 
graphene, with its gapless chiral linear 2D electron-hole 
Dirac band dispersion, the resistivity (conductivity) will 
increase (decrease) monotonically with decreasing car- 
rier density near the charge neutrality (Dirac) point pro- 
vided the system is sufficiently pure, i.e., ordered, with 
very little residual background charged impurity disor- 
der. Not only will the Dirac point resistivity be anoma- 
lously large in high-purity graphene, the transport behav- 
ior itself will be insulating-like at the charge neutrality 
point with the resistivity increasing monotonically with 
decreasing temperature! On the other hand, as the car- 
rier density increases, the resistivity will decrease with 
the eventual restoration of the metallic behavior mani- 
festing a weakly temperature dependent resistivity above 
a non-universal crossover density which would depend on 
the residual background charged impurity disorder and 
the temperature. At high density, far away from the 
charge neutrality point with vanishing average charge 
density, the high-purity graphene sample would behave 
in a perfectly normal manner manifesting very high mo- 
bility (and very long mean free path) consistent with 
the highly ordered nature of the system with very lit- 
tle residual Coulomb impurity scattering. We dub this 
strange dichotomy where decreasing disorder drives the 
graphene layer into an effective insulating state at low 
carrier density near the charge neutrality point, while 
maintaining very high mobility at high carrier density 
consistent with its low disorder, the phenomenon of "dis- 



order by order" . We emphasize that our predicted dis- 
order by order phenomenon is not a T = quantum 
phase transition as in an Anderson or Mott transition, it 
is a transport crossover phenomenon manifesting itself as 
an effective density-tuned metal-insulator transition. In 
particular, quantum localization plays no role in our the- 
ory which is developed entirely within the semiclassical 
Boltzmann transport model neglecting all quantum inter- 
ference corrections. The disorder by order phenomenon 
arises from an interplay among charged impurity disor- 
der, density inhomogeneity (the so-called "electron-hole" 
puddles^Ti^) , and the peculiar gapless linear chiral band 
dispersion of graphene. Our predicted novel semiclassi- 
cal phenomenon would dominate low-density transport 
in ultrapure graphene samples as long as quantum inter- 
ference induced localization corrections are small, i.e., in 
the effective high-temperature semiclassical regime where 
the inelastic phase breaking length is comparable to or 
smaller than the elastic transport mean free path. 



We believe that our predicted graphene "disorder by 
order" phenomenon has recently been experimentally ob- 
served by Ponomarenko et alii-, who, however, interpret 
their observation as the manifestation of a density-tuned 
metal-insulator localization transition. All aspects of 
the experimental data reported in Ref. [l|, in particu- 
lar the density and the temperature dependence of the 
measured conductivity, agree spectacularly well with our 
predictions, and we therefore contend that the observa- 
tion in Ref. [l| is a direct experimental verification of our 
predicted "disorder by order" phenomenon. Particularly 
germane in this context is the fact that the transport data 
of Ref. [l| were taken in the relatively high-temperature 
(10 — 100 K) regime where quantum interference effects 
are not important, and our semiclassical transport theory 
shoidd apply [See Appendix [B] for details]. 

The paper is organized as follows. In Sec. [TTl the basic 
transport theory and the numerical results are presented. 
In Section [ml we discuss the results compared to exper- 
imental data, and we conclude in Sec. IIVI 
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II. THEORY AND NUMERICAL RESULTS 

We first provide a simple physieal pieture underly- 
ing the "disorder by order" phenomenon. Let us as- 
sume that the graphene sample is pristine with essen- 
tially no baekground random charged impurities so that 
the effective transport relaxation or scattering time r is 
very long, leading to very high (low) conductivity (re- 
sistivity) at an electron density of n. The conductiv- 
ity is given (at T = 0) within the Boltzmann trans- 
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port theorji^i^ by cr = ^-j^D{Ep)t{Ef), where vp is 
the constant graphene Fermi velocity defining its linear 
band dispersion, D{Ep) cx ^fnjvp is the graphene den- 
sity of states at the Fermi energy Ep (x vp^/n. The 
scattering time t{Ep) for the screened Coulomb scatter- 
ing due to random background charged impurities has 
been calculated^"— in the literature, giving r oc ^/n, 
which leads to the now- well-known formulai^ for graphene 
conductivity due to random charged impurity scattering 
given by 



where n,; is the effective baekground 2D concentration 
of the random charged impurities (including its location 
and strength) whereas A is a constant which depends on 
the dielectric environment of the system (e.g., substrate)- 
for graphene on Si02 (h-BN), A » 20 (26). The resuh 
given in Eq. [1] and the underlying general theory for 
graphene carrier transport have been well-verified exper- 
imentally in the literature^. An immediate consequence 
of Eq. [U interpreted naively, is that the graphene resis- 
tivity p(= 1/cr) diverges as at the Dirac point where 
the carrier density vanishes by virtue of the vanishing 
density of states at the Dirac point. 

In reality, however, this divergent Dirac point resistiv- 
ity (or equivalently, vanishing conductivity) is not ob- 
served experimentally in real graphene samples, which 
manifest a conductivity saturation phenomenon at low 
carrier density (|n| < n*—) with an approximate mini- 
mum conductivity plateau amm An* /ui, where n* is 
a characteristic disorder-dependent density^i^. This con- 
ductivity minimum phenomenon was already apparent in 
the pioneering graphene experiments by Novoselov and 
Geiro^, and was later studied extensively quantitatively^, 
and is now accepted as arising from the charged impu- 
rity disorder induced inhomogeneous electron-hole den- 
sity puddles which dominate the graphene landscape^ at 
low carrier density. These puddles of strong real-space 
density inhomogeneities arise from the low-density fail- 
ure of screening of the individual charged impurities with 
electrons/holes preferably accumulating near/far from 
individual discrete impurities depending on the sign of 
the impurity charge^^. Thus, as the gate voltage de- 
creases, the average density decreases, but electron-hole 
puddle formation leads to an effective saturation of the 
conductivity at some low sample-dependent minimum 



value. The inhomogeneous puddles simply cut off the 
p ^ 1/n behavior of graphene resistivity for n "^n* since 
the real 2D density across the graphene sample never 
vanishes although the average density does, allowing for 
percolating transport through the electron-hole puddles 
at the charge neutrality point^. 

What would happen if the electron-hole puddles are 
somehow eliminated or suppressed in the system? Within 
the semiclassical Boltzmann picture, the resistivity will 
become very large as the average density is decreased by 
lowering the gate voltage since the puddles leading to 
the low-density conductivity saturation phenomenon no 
longer exist! This is a direct (and dramatic) manifesta- 
tion of the gaplessness of graphene, and cannot happen 
in the semiconductors with band gaps. 

The easiest way to eliminate (or suppress) the puddles 
is, of course, to reduce the environmental charged impu- 
rity density {rii) which induces the puddles to start with. 
But such a low-disorder system will necessarily manifest 
very low resistivity (since p (x Ui) at high carrier den- 
sity (p cx 1/?T.), but very high resistivity near the Dirac 
point since n — > 0. If the puddles disappear completely, 
the resistivity will diverge as 1/n as the carrier density 
decreases. Therefore, the disorder-by-order phenomenon 
is peculiar to gapless graphene with its linear dispersion. 
It is obvious from the above physically-motivated discus- 
sion based on a qualitative extension of existing results in 
the literature^iii^iiS that this counter-intuitive "disorder- 
by-order" phenomenon would be more apparent if the 
inhomogeneous electron-hole puddles can be further sup- 
pressed around the Dirac point by applying an external 
screening potential through a gate which would screen 
out the puddles, as has been successfully done in Ref. 

The above-discussed semiclassical "disorder-by-order" 
phenomenon has recently been observed in the experi- 
ment of Ponomarenko et alX, who reported monotonie 
increase of the graphene resistivity with decreasing car- 
rier density in an ultrapure sample on h-BN substrate. 
This remarkable resistivity enhancement with decreasing 
density occurs only in the presence of a second nearby 
high-density graphene layer which screens out the pud- 
dles, thus avoiding the "minimum-conductivity" satura- 
tion phenomenon around the Dirac point. 

In Fig. [I] we present our theoretically calculated trans- 
port results as a function of average carrier density for 
the experimental situation studied by Ponomarenko et 
ali^. The different colors in Fig. [1] correspond to dif- 
ferent temperatures whereas different panels correspond 
to different electron-hole puddle configurations charac- 
terized by the disorder induced potential fluctuation 
parameter 's' where 's' corresponds to the root mean 
square potential fluctuations in the probability distribu- 
tion function P{V) for the impurity-induced disorder, 
P{V) ~ e"^ /V2Trs^, assumed to be Gaussian for 
simplicity - the Gaussian approximation is very accurate 
compared with the realistic numerical calculation a^'^*^ of 
P(V). The inhomogeneity parameter 's', which depends 
on the impurity disorder in the system, is directly con- 
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FIG. 1. Calculated p(n) at different temperatures T for 
rii = 16 X 10^" cm~^. (a) The potential fluctuation parame- 
ter s = 22 meV. (b) s = 10 meV. (c) s = 6.0 meV. (d) s = 1.0 
meV. fe) s = 0.1 meV. (f) Experimental data as shown in 
Ref. [l[. Comparison between (b) and (f) indicates a poten- 
tial fluctuation parameter s ~ 10 mev in the experimental 
sample^. 



nected to the room-mcan-squarc density fluctuation n„ns 
in the inhomogeneous electron- hole puddles. The precise 
relationship between s and Urms can only be obtained 
through a full numerical self-consistent calculation^ii^, 
but within a simple mean- field theory rirms oc s^. We 
note that rirms ~ n* defines the cut-off for the minimum 



s around the Dirac point 
0, PCNP = l/o-„„;„ diverges 



conductivity amin ^ »^rm; 

as discussed above. As s 
as in the mean field theory. We note that the poten- 
tial fluctuation s (or equivalently, the root mean square 
density fluctuation in the puddles) is being controlled 
by external gating through the second graphene layer in 
Ref. and the new feature of Ref. not achieved 
before, is that s could be made very small. 

For the sake of comparison, we have reproduced in 
Fig. I^f) the corresponding experimental results from 
Ref. The agreement between our calculated theoret- 
ical results and the experimental data is striking: In the 
presence of substantial (vanishing) electron-hole puddles 
characterized by larger (smaller) values of the disorder 
fluctuation parameter s, the calculated p{n) saturates 
(increases monotonically) at lower average carrier den- 
sity exactly as observed experimentally in the presence 
(absence) of the puddles. We emphasize that the the- 
oretical results obtained in Fig. [T] use exactly the same 
parameters for all cases except for varying the value of the 
potential fluctuation parameter 's' as given in the figure 
(which mimics the effect of the suppression of s through 
screening by the second graphene layer in Ref. [l|). One 
interesting prediction of our theory is that p{n) never 
truly diverges as 1/n in our theory as n — > (unless 



s = exactly, which is unphysical) since there is always 
a low-density cut-off n* defining the conductivity mini- 
mum regime with n* decreasing with decreasing s. This 
low-density cut-off (and the corresponding "maximum re- 
sistivity" Pmax) depends strongly on the puddle param- 
eter 's' - the suppression of 's' dramatically increases 
(decreases) Pmaxin*). 

Before discussing our results for the temperature de- 
pendence of the resistivity p{T), we briefly discuss our 
transport theory^ili^iiSiiii for graphene conductivity in 
the presence of electron- hole puddle induced strong den- 
sity inhomogeneity. The conductivity is obtained by us- 
ing the effective medium theory (EMT) by solving the 
integral equation 



cr(n) - (TEMT pr 1 „ 

an—— F[n\ = 

a(n) + UEMT 



(2) 



where uemt is the effective conductivity of the sample 
and cr(n) is the density n{r) dependent local conduc- 
tivity with the carrier density n having the distribution 
Pin) = exp[— (n — no)^/2ri^„g]/-\/27rn2^^ defining the 
electron-hole puddles - here tiq is the average density 
defined by the external gate voltage (i.e., no = at 
the charge neutral Dirac point) and n™is is the root- 
mean-square density fluctuation due to the existence of 
density inhomogeneity associated with the puddles. We 
calculate o'(n) using the finite-temperature Boltzmann- 
RPA transport theory using screened random quenched 
charged impurity centers (of 2D concentration ni) in 
the environment as the resistive scattering mechanism. 
The Boltzmann transport theory, which has been de- 
scribed in details elsewhere^iiiii^ , includes five distinct 
temperature-dependent contributions: (1) thermal acti- 
vation of electron-hole occupancy (i.e. thermal excita- 
tion of electrons from the valence band to the conduction 
band); (2) finite temperature thermal averaging around 
the Fermi surface according to the Fermi distribution 
function; (3) the thermal activation of carriers over the 
potential fiuctuations associated with the electron-hole 
puddles; (4) finite temperature screening by the carriers 
themselves; (5) phonon effects (which are straightforward 
to includei^, but are neglected here since electron-phonon 
coupling is weak in graphene). The thermal effects (1)- 
(3) above produce 'insulating' temperature dependence, 
i.e., the temperature-dependent resistivity p{T) increases 
with decreasing T, whereas the last two effects lead to a 
'metallic' p{T) decreasing with decreasing temperature. 
All the thermal effects are suppressed with increasing car- 
rier density (or more precisely, increasing Ep), and they 
arc the strongest at the charge neutral Dirac point (where 
the nominal Ep vanishes). We note that in our figures, 
p = <^EMT whereas the density n = n^, i.e., the average 
density. 

In Fig. [2] we depict our calculated p{T) at fixed n in- 
cluding puddle effects (characterized by the parameter 
s) within the Boltzmann transport theory as described 
above (and clsewher o^^'^^ ). The calculated temperature 



4 




10 100 10 100 




.5 



T(K) n(10"cm^2) 

FIG. 2. Calculated p{T) for different values of the potential 
fluctuation parameter s. (a) and (b) correspond to n = 0, 
i.e., the Dirac point, (c) and (d) correspond to high carrier 
density n = 5x10^" cm~^. (a) and (c) rij = 9xI0^° cm"^ with 
mobility at high carrier density fj, = 10^ cm^/Vs. (b) and (d) 
rii = 16 X 10^'' cm~^ with mobility at high carrier density fi — 
6 X 10* cm^/Vs. The decrease of p for T > 100 K in (c) and 
(d) is a real effect arising from the Fermi surface averagingi^ 
which becomes quantitatively important at these densities - 
this effect is suppressed by phonon scattering which comes 
into play for T > 100 K. (e) for different values of T are 
shown for s — 5 meV and = 9 x 10^" cm^'^. 



dependent resistivity is identical to the experimental ob- 
servations of Ref. [l|, as reproduced in our Fig.[IJf), with 
the Dirac point resistivity increasing strongly with low- 
ering temperature for smaller values of s (i.e., when the 
potential fluctuations associated with the electron-hole 
puddles are strongly suppressed) whereas the resistivity 
at high carrier density manifesting almost temperature- 
independent behavior. We emphasize, however, that 
within our semiclassical transport theory, in spite of the 
very strong increase of the Dirac point p{T) with lower- 
ing T for small s, eventually p{T) saturates at some large 
s-dependent (and ni-dependent) value at low enough 
T <C s even at the Dirac point - empirically we find that 
p{T) at the Dirac point saturates for T < s/5. Then, 
the Dirac point behavior of p{T) for different values of s 
is qualitatively similar if plotted as a function of ksT/s. 
Therefore, one clear prediction of our semiclassical theory 
is that at low enough temperatures the experimentally 
measured p{T) will always saturate even at the Dirac 
point, but the crossover temperature for this saturation 
would be very low when the electron-hole puddles are 
strongly suppressed. Thus, p{n, T) in our theory would 
behave very much like a density-tuned metal-insulator 
transition, as observed in Ref. for all practical pur- 
poses except that the temperature dependence would be 



a power law (and not an exponential as in a strongly 
localized system), exactly as seen in Ref. 



III. DISCUSSION 

One point we should make clear here about the "disor- 
der by order" phenomenon is that it arises from the sup- 
pression of the electron-hole puddles associated with the 
density inhomogeneity around the Dirac point induced 
by Coulomb disorder, and not simply from an impurity- 
induced collisional broadening effect in the graphene den- 
sity of states. In particular, one may wonder whether 
the impurity broadening effect leading to a finite carrier 
density of states at the Dirac point could by itself lead 
to the disorder by order phenomenon since this would 
imply that a cleaner system would have a lower density 
of states at the Dirac point and hence a higher resistiv- 
ity. While this is certainly true in general, i.e. when 
a material exhibits a density of states minimum, addi- 
tion of disorder leads to an increased density of states 
near this minimum due to smearing which then leads to a 
higher Drude conductivity, this is not the primary opera- 
tional mechanism for the disorder by order phenomenon. 
As we show in Appendix |X] of this paper, such an en- 
hanced disorder-broadened density of states indeed leads 
to higher conductivity near the Dirac point for more dis- 
ordered graphene samples, but the conductivity always 
vanishes at the Dirac point {n = 0) in this situation as 
long as Coulomb disorder (i.e. random charged impuri- 
ties in the environment) is present in the system inde- 
pendent of whether the density of states is zero or finite 
at the Dirac point. In our proposed order by disorder 
mechanism, by contrast, the conductivity is always fi- 
nite at the Dirac point because of the disorder-induced 
inhomogeneous puddles except that this finite minimum 
Dirac point conductivity arising from the puddles is much 
lower for the more ordered systems than it is for the more 
disordered systems while the reverse is true at high den- 
sity. Thus, the presence (suppression) of puddles is the 
important physics in our disorder by order mechanism, 
not the disorder-induced smearing of the graphene den- 
sity of states around the Dirac point. We discuss the issue 
of the density of states smearing effect on the graphene 
conductivity in Appendix 1X1 

Before concluding, we now provide a critical discus- 
sion of our theory as applied to the experimental obser- 
vations of Ref. m. This issue is discussed in details in 
Appendix [B] First, our results provide an excellent de- 
scription of the experimental observations with respect 
to the dependence of the measured resistivity as a func- 
tion of carrier density, temperature, and disorder associ- 
ated with the electron-hole puddles. Indeed, the agree- 
ment between our theory and the experimental data is 
striking, lending credence to our claim that the obser- 
vation of Ref. [l[ is an experimental verification of our 
predicted "disorder by order" phenomenon in monolayer 
graphenoii. In Ref. the potential fluctuations asso- 
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ciated with electron-hole puddles in the graphene layer 
were suppressed by a second close-by graphene layer with 
a very high carrier concentration. We have verified by 
direct numerical simulations that such a suppression can 
indeed be caused by the second high-density graphene 
layer separated by a distance d acting as a gate which 
screens the potential fluctuations reducing them substan- 
tially below their pristine value arising from the random 
charged impurity distribution. A simple electrostatic an- 
alytic calculation shows that the suppression of s would 
be approximately by a factor l/{kpd) when kpd <C 1. 
For n ^ 10^° cm~^ and d ~ 1 nm, s could thus be sup- 
pressed by a factor as large as 50! Thus our basic picture 
of the suppression of the potential fluctuation parame- 
ter 's' leading to the disorder by order phenomenon at 
the Dirac point is applicable to Rcf. In Rcf. it 
is also found that the disordcr-by-order phenomenon is 
itself suppressed by the application of a weak magnetic 
field B ^ 10 mT. Although the inclusion of a magnetic 
field in the theory is beyond the scope of the current 
work, we mention that a magnetic field B ~ 10 mT cor- 
responds to a minimum Landau level energy separation 
~ 4 meV which is comparable to Ep for n ^ 10^ cm~^. 
Thus, a 10 mT magnetic field is not a weak field at the 
Dirac point, which would nonperturbatively modify the 
physics, considerably suppressing the disorder by order 
phenomenon. 

We believe that the disorder by order phenomenon 
should occur in any Dirac material with chiral gapless lin- 
ear energy spectrum as long as quantum interference ef- 
fects are negligible, and as such, we predict the existence 
of the same phenomenon in the 2D surface transport^^ 
in 3D topological insulators provided that the surface 
puddles are suppressed in the system and the temper- 
ature is not too low. In fact, all gapless semiconduc- 
tors will manifest the disorder by order phenomenon 
if impurity-induced potential fluctuations can be sup- 
pressed in the low carrier density regime. Indeed, we be- 
lieve that some earlier graphene experiments^^ observing 
anomalous temperature dependence of the Dirac point re- 
sistivity are observing exactly the same disordcr-by-order 
phenomenon as reported in Ref. [l[ except that the au- 
thors of Ref. [l^ interpret their observations as ballistic 
transport whereas Ponomarenko et alX invoke Anderson 
localization! 



IV. CONCLUSION 

We conclude by emphasizing that the recently ob- 
served intriguing phenomenoii^ of monotonically increas- 
ing graphene resistivity with decreasing carrier density in 
graphene samples of very high purity (with very high 
mobilities at high carrier densities) most likely arises 
not from Anderson localization, but from the semiclas- 
sical "disorder by order" phenomenon proposed in the 
current manuscript. This phenomenon arises from the 
suppression of electron-hole puddles in the system by a 



near-by screening layer which then induces the system to 
show its intrinsic Drude behavior of the resistivity be- 
ing inversely proportional to the carrier density down to 
much lower carrier densities without being cut off by the 
puddle-induced (and well-known2iiii£iii) "graphene con- 
ductivity minimum" mechanism. The qualitative differ- 
ence between our "disorder by order" mechanism and 
localization is that in our case, the conductivity is al- 
ways finite, eventually being cut off by remnant puddles 
in the system at much lower carrier densities determined 
by the details of disorder and screening by the second 
layer whereas in for localization the conductivity is truly 
zero at T = 0. The other qualitative difference is that 
the predicted temperature dependence of the resistivity 
at a fixed low carrier density near the Dirac point in our 
mechanism is a power law whereas it must be exponen- 
tial in the case of Anderson localization. The observed 
temperature and density dependence of the low-density 
resistivity in Ref. [l| is consistent with our predictions, 
and hence we believe that Ref. [l[ is manifesting the "dis- 
order by order" phenomenon, not Anderson localization. 
When and how graphene can manifest localization (as op- 
posed to antilocalization, which is the expected behavior 
for graphene) is an interesting question for the future, but 
we believe that the experimental temperature range must 
be much lower than that used (20 — 100 K) in Ref. [ij to 
see any quantum interference induced localization effects 
since the inelastic phase coherence length is simply too 
short (Appendix |B|) for localization effects to manifest at 
higher temperatures. We urge transport experiments in 
high quality graphene at very low temperatures (< 1 K) 
to discern localization/antilocalization versus scmiclassi- 
cal "disorder by order" puddle effects. 



ACKNOWLEDGMENTS 

This work is supported by US-ONR-MURI and NRI- 
SWAN. One of the authors (SDS) gratefully thanks An- 
dre Geim for stimulating discussions and helpful com- 
ments. 



Appendix A 

The density of states (DOS) of disordered graphene is 
calculated within the self-consistent Born approximation 
(SCBA). In calculating the DOS only a short range dis- 
order potential is considered. The T = conductivity is 
calculated with a simple formula 

a{n) = ^D{Ep)r{Ep), (Al) 

where D{Ep) is the DOS at Fermi level, and t{Ep) is 
the transport scattering time. Note that the scattering 
time is calculated with the DOS of the bare band and 
with two different disorders, short range potential and 
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long range Coulomb potential. Thus, the scattering time 
is given by 



1 



1 

To 



(A2) 



where tc (tq) is the scattering time due to the long-range 
Coulomb impurities (short range impurities) and they are 
given by in the Boltzmann transport approach 

— Db{EF)\Vc{kF)\^ oc Ep\ 

- cx Dt,{EF)\VQ{kF)\'' cx Ef, (A3) 

To 

where Db{E) oc i? is the DOS of bare band. In low 
density limit, Ef ^ 0, the scattering rate by Coulomb 
impurities dominates and the total scattering time be- 



cx n 



1/2. Thus, 



haves as t{Ef) oc Ef, i.e., t{Ef) 
though the DOS is finite at a; = the conductivity, 
a{n) oc D{Ef)t{Ef) oc y/n. as n ^ 0. 

We calculate the disorder-broadened DOS for graphene 
following Refs. and then calculate the conductiv- 

ity following Eqs. (jAip - (jA3p above. Our results for the 
broadened DOS and the resulting conductivity are shown 
in Figs. 13] and m The important points to note are :(1) 
the disorder-broadened DOS leads to an enhanced (sup- 
pressed) conductivity at low (high) carrier densities for 
more (less) disordered systems in agreement with Ref. [l| ; 
(2) but the Dirac point conductivity is always zero in- 
dependent of whether the DOS is smeared by disorder 
or not and thus the smearing of the DOS by itself can- 
not be the explanation for the existence of the minimum 
conductivity plateau in graphene which necessitates the 
existence of electron-hole puddles in the system. 

We can consider the DOS of disordered graphene 
(rather than the DOS of the bare graphene as done 
above) to calculate the scattering time. Now we have 
a finite DOS as — >■ 0. Then with the same approach as 
Eq. (|A3I) we have the scattering times as Ef — >■ 



— oc D(EF)\Vc(kF)\^ oc E-^ 
Tc 

— oc D{EF)\Vo{kF)\^ oc E"p. 

To 



(A4) 



Similarly we have r oc Ep oc n, and cr{n) oc n as n — >■ 0. 

Thus, the DOS smearing by disorder always produces 
zero conductivity at the Dirac point (n = 0) in graphene 
although close to the Dirac point the smearing of the 
DOS does indeed lead to an enhanced conductivity as 
shown in the figures. Although our results are shown 
within the SCBA theory of the DOS smearing, the qual- 
itative findings are the same within the simpler Born ap- 
proximation where the disorder broadening of the DOS 
does not lead to a finite DOS at zero density. Thus, 
the DOS smearing by disorder cannot be an explanation 
for the graphene finite minimum conductivity around the 
Dirac point which arises from the Coulomb disorder in- 



^0.6 
Q 

'3 0-4 

Q 

0.2 



n'y„2=x(eVA)2 


/ 

/ 


-(a) 






50 






r , \ , \ , 


1 



0.2 
^0.15 

Q 

3~ 0.1 
Q 

0.05 



n,Vo^=x(eVA)2^ 




;(b) ^ 




















1,1,1 


, 1 , 



0.2 



0.4 
CO/Ec 



0.6 0.8 



0.02 0.04 0.06 0.08 

(0/Ec 



0.1 



FIG. 3. Calculated the density of states within SCBA for 
different strengths of (5-range disorder potential Vb, ndVo = 0, 
10, 50, 100 (eVA)^ where Ud is the short range impu- 
rity density. Here Ec ~ hvpkc is the cut-off energy with 
kc ~ 1/a where a = 1.42A is the lattice constant. Do = 
gsgvEc/{2Tv{hvF)^)- (b) shows the same results of (a) at low 
energy regime. The DOS at a; = is given by the formula, 
D{0)/Do = \nEo/{2nEo), where Eq = E^j^e^-^l^i - 1 and 
7 = ndK?/(2(ft«F)'). 
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FIG. 4. (color online) Calculated conductivity as a func- 
tion carrier density. Parameters of graphene on h-BN are 
used. The conductivity is calculated with Eq. lAll a = 

^-^D(Ef)t{Ef), where D{Ef) is the DOS at Fermi level 
which is given in Fig. [3l and t{Ef) is the transport scatter- 
ing time calculated with DOS for the bare band. The con- 
ductivity is calculated for a fixed long range charged impurity 
Ui = 10^1 cm~^ , but for different short range disorder poten- 
tials as given in the figure. Both the DOS and conductivity 
are calculated with the same short range disorders. Note that 
the conductivity approaches zero as the carrier density goes 
to zero, (b) shows the resistivity p — 1/a calculated with the 
same parameters of (a) as a function carrier density in log-log 
scale. 



duced density inhomogeneity and electron-hole puddles 
in the system. 



Appendix B 

In this Appendix, we tackle two important (and inter- 
related) questions: (i) What about Anderson localiza- 
tion? (ii) Can our Boltzmann theory be valid for p > 
h/e^l We discuss these questions in some details below. 

First, in the absence of inter- valley scattering graphene 
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should manifest anti-localizationi^ behavior for which 
there is some experimental evidence in some situationsi^. 
The presence of inter- valley scattering would restore the 
usual localization behavior—. But the crossover from 
the anti-localization behavior to the localization behavior 
should occur at an extremely low temperature compara- 
ble to the inter-vallcy scattering rate which is likely to be 
much less than 1 K. It is therefore reasonable to assume 
that localization effects are operational^^ at much lower 
temperature scales than the results presented in this pa- 
per and that used in the experiments of Ref. [l| . In addi- 
tion, the measurements of Ref. [ij are performed at high 
temperatures (10— 100 K) where quantum interference ef- 
fects are strongly suppressed because the inelastic phase 
breaking length is shorter than the elastic mean free path. 
The validity of the Boltzmann theory is, as a matter of 
principle, completely independent of the resistivity of the 
system as long as quantum localization/interference ef- 
fects can be neglected. Thus, our semiclassical theory 
does not simply care about the condition p > h/e^ as 
long as quantum interference corrections are small which 
is true at high temperatures. 

We discuss below the extent to which our semiclassi- 
cal disorder-by-order transport phenomenon applies to 
the recent experimental observations of Ponomarenko et 
alr^, who have interpreted their measured graphene resis- 
tivity p{n, T) in terms of an Anderson localization driven 
metal-insulator transition at low carrier density. 

We believe that the data of Ref. [l[ are inconsistent 
with quantum interference induced localization transi- 
tion for the following reasons: (1) the temperature regime 
explored ini is sufficiently high (10 — 100 K) where local- 
ization effects should be strongly suppressed since the in- 
elastic phase breaking length is comparable to or shorter 
than the elastic mean free path making quantum inter- 
ference irrelevant; (2) the temperature dependence of 
the measured resistivity in Ref. [H is not exponential 
and is thus inconsistent with Anderson localized insu- 
lating behavior; (3) the density dependence of the mea- 
sured resistivity in Ref. [ij follows our predicted semi- 
classical 1/n behavior at low density (before it satu- 
rates at a very low temperature-dependent cut-off den- 
sity near the Dirac point exactly as our semiclassical 
theory predicts)- this is hard to reconcile with Ander- 
son localization; (4) there is no observable weak local- 
ization (or anti-localization) behavior manifesting in the 
data at high density where the resistivity is low (and 
weak localization corrections should be discernible if in- 
deed a density-tuned metal-insulator strong localization 
phenomenon is taking place); (5) graphene can manifest 
true localization (even at T = 0) only in the presence 
of intervalley scattering which is generally known to be 
extremely weak (in the absence of intervalley scattering 
and trigonal warping, graphene can only manifest anti- 
localization- the experimental data of Ref. [l[ do not di- 
rectly exhibit any signature of strong inter-valley scat- 
tering, for example, the high-density resistivity is com- 
pletely consistent with the intervalley scattering being 



weak); (6) the temperature scale for the strong local- 
ization transition should be lower than the characteris- 
tic (very low) energy scale for the intervalley scattering, 
but the experimental data of Ref. [l| manifest insulating- 
like behavior already at rather high (^ 10 K) tempera- 
tures, casting doubts on the whole metal-insulator tran- 
sition picture; and finally (7) as emphasized in our main 
manuscript, our semiclassical theory is in excellent agree- 
ment with the observed experimental data of Ref. [l[ , and 
the observed data manifest no direct signature of quan- 
tum interference effects. 

Having argued above that the experimental observa- 
tions of Ponomarenko et alX arc unlikely to be aris- 
ing from a quantum strong localization induced metal- 
insulator transition, we now give the reasons for our 
firm belief that Ref. [l[ is a direct manifestation of our 
predicted semiclassical "disorder by order" phenomenon: 
(1) our theoretical results provide an excellent quantita- 
tive description of the experimentally measured resistiv- 
ity; (2) the density and the temperature dependence of 
the resistivity in Ref. [l| follows precisely the asymptotic 
p ^ 1/n and p ^ 1/T^ behavior predicted in our the- 
ory; (3) we predict the correct disorder dependence, i.e., 
lower the puddle induced inhomogeneity or potential fluc- 
tuation, stronger is the low-density insulating behavior; 
(4) our theory provides a clear qualitative and quanti- 
tative explanation for why the inhomogeneous electron- 
hole puddles must be suppressed in order to observe the 
apparent low-density insulating behavior; (5) our the- 
ory can explain the data both in the presence and in 
the absence of the puddles (with the low-density insu- 
lating behavior being respectively absent and present) 
whereas in the quantum localization scenario no explana- 
tion is available for the generic situation in the presence 
of puddles- one must assert ad hoc that the presence of 
puddles somehow hinders the localization effect; (6) in 
the high temperature regime explored in Ref. [ij, quan- 
tum interference effects are strongly suppressed making 
our semiclassical theory applicable - we emphasize that 
quantum interference cannot occur if the phase coher- 
ence length is short as it is in the temperature regime of 
Ref. g. 

We come to the purely theoretical (rather than em- 
pirical) question of the applicability of our semiclassical 
theory to the experimental situation of Ref. [l[ where the 
measured resistivity is high [p > h/e^) so that the well- 
known localization condition kpl^ > 1 (which is equiv- 
alent to the p > h/e^ condition written out in terms of 
Fermi wavcvector and mean free path) is violated, and 
the low-density insulating phase in Ref. [l| actually cor- 
responds to kple < 1 regime, where Ig is the elastic (or 
transport) mean free path, given by 1^ = vpT, where 
T is the appropriate finite-temperature elastic scatter- 
ing time due to Coulomb disorder. It is indeed true that 
kple ^ 1 is approximately the condition ("the loffe-Regel 
criterion" ) for localization effects to become important 
in an electronic system with the regime kplg < 1 be- 
ing dominated by quantum interference induced localiza- 



8 



tion effects where semiclassical transport theory should 
be nominally invalid since localization corrections to the 
resistivity should become larger that the semiclassical re- 
sistivity itself. This is, however, true only at T = (or at 
extremely low temperatures) where k 3> le with li being 
the inelastic phase coherence length for the electrons (at 
T = 0, — oo, and the condition li ^ Ig is trivially 
satisfied). Thus, quantum localization requires a neces- 
sary condition {li ^ l^) so that quantum interference 
is operational and a sufficient condition (If > 1^, where 
If ^ 1/kp is the Fermi wavelength). Without the neces- 
sary condition {li 3> h) being satisfied, quantum inter- 
ference is simply not operational in a given experimental 
situation. (This is the reason that quantum localization 
phenomenon is typically experimentally studied at very 
low mK range of temperatures, not at T > 10 K.) 

The necessary condition [li 3> le) for quantum interfer- 
ence (i.e., localization) to be operational is not satisfied 
by the experimental conditions (T = 10 — 100 K) used in 
Ref. [l|. In particular, we have calculated^! the inelas- 
tic phase breaking length li due to electron-electron and 
electron-phonon interactions in graphene, finding that 
typically k < le in the higher temperature (> 20 K) 
regime used in Ref. [l[, and even at the lowest measure- 
ment temperature in Ref. [l|, T ~ 10 K, li < le- Our the- 
oretical results for li will be published elsewhere^i, but 
we mention that our finding that li < le in the experi- 
mental regime of Ref. [l[ is consistent with the available 
direct experimental measurements of h in graphene. For 
example, Ref. [l| specifically quotes li ^ /im at liquid he- 
lium temperature (T ~ 4 K). which translates into li ^ 
10 - 400 A for T = 100 - 20 K since k - T'"^ in graphene 
due to the dominant electron-electron interactional. The 
actual values of U may even be smaller because of dis- 
order effects and electron-phonon interactions. Direct 
measurements of h in graphene^^ are consistent with our 
estimate that h < le in most of the temperature regime 
explored in Ref. [l|. Thus, quantum interference phe- 
nomena are unlikely to be playing any role in Ref. [l| 
except perhaps at the lowest temperatures (T '--^ 10 K). 

It will be desirable for the experimental regime of 
Ref. [H to be pushed down well below 10 K (T < 1 
K) where the necessary condition for quantum interfer- 
ence to be operational, namely, li ^ le-, should apply. 
We expect our theory to fail in this low-temperature 
interference-dominated puddle-free regime, and it will be 
extremely interesting to study the deviations of the ex- 
perimental data from our semi-classical theory in this 
low temperature (T < 1 K) regime so as to learn about 
the nature of graphene localization. We believe, for rea- 
sons discussed above, that the current experimental data 



of Ref. [l| fall in the intriguing regime of k < le and 
kple ^ 1 where our semiclassical theory still remains 
valid. 

We point out a technical issue which makes our the- 
ory even more applicable to the measurements in Ref. [l[ 
than the above discussion implies. The elastic scatter- 
ing length le that li should be compared to (so as to 
check if li ^ le condition is satisfied for quantum in- 
terference to be relevant) is not the net mean free path 
defining the resistivity, but only the inter- valley scatter- 
ing length liy le- This is because most of the elas- 
tic resistive scattering in graphene is intravalley scatter- 
ing which cannot lead to Anderson localization (in fact, 
it should induce anti-localizationi^) . Since the actual 
inter- valley scattering is very weak, according to Ref. [l[ , 
liv ~ 0.1 — O.Sfiui, we conclude that the experimental 
situation in Ref. [l[ corresponds entirely to the regime 
k < ky in the T = 10 - 100 K range since k T"^ 
and therefore k < 1000 A throughout the temperature 
range studied in Ref. [l| . Basically, the experimental con- 
ditions in Ref. [l| correspond more to the infinite tem- 
perature limit than the zero temperature limit as far as 
quantum interference effects are concerned and as such 
our semiclassical theory should apply well to Ref. [l[ . We 
therefore believe that the current data of Ref. [l| are well- 
explained by our theory to be a semiclassical disorder by 
order phenomenon. 

We conclude by noting that in addition to the semi- 
classical disorder by order and the localization induced 
metal insulator transition phenomena, in principle there 
is a third possible explanation for the low-density insu- 
lating behavior observed by Ponomarenko et alX. This is 
the possibility of a spontaneous interaction-induced gap 
formation at the Dirac point^, which would of course 
lead to the insulating behavior when the chemical poten- 
tial approaches the gap at low density (similar to ordi- 
nary semiconductors). This gap scenario is unlikely to 
be operational in the experiment of Ref. [l[ because (1) 
typically the predicted spontaneous gap is very small, 
much less than the temperature range (T = 10 — 100 
K) explored in Ref. [ij, and (2) the temperature depen- 
dence of the resistivity in Ref. [l| is a power law (roughly 
T~^) which is inconsistent with the exponentially acti- 
vated resistivity expected for a gapped system. Again, 
experiments need to be carried out at much lower tem- 
peratures (below 1 K) for a thorough investigation of 
this issue of the spontaneous gap formation, which has 
recently been controversial in the literature with different 
theories claiming both the existence and the nonexistence 
of an interaction-induced graphene gap formatioiiSl. 
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